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ABSTRACT 

A  complete  asymptotic  development  of  the  Stirling  numbers 
S(N,K)   of  the  second  kind  is  obtained  by  the  saddle  point  method 
previously  employed  by  Moser  and  Wyman  [Trans.  Roy.  Soc.  Cannad., 
49(1955)  49-54]  and  deBruijn  [Asymptotic  Methods  in  Analysis, 
North-Holland  Publishing  Co.,  Amsterdam,  (1958)  102-109]  for  the 
asymptotic  representation  of  the  related  Bell  numbers. 


ASYMPTOTICS  OF  STIRLING  NUMBERS  OF  THE  SECOND  KIND 


W.  E.  Bleick  and  Peter  C.  C.  Wang* 

Mathematics  Department 
Naval  Postgraduate  School,  Monterey,  California 


A  complete  asymptotic  development  of  the  Stirling  numbers 
S(N,K)   of  the  second  kind  is  obtained  by  the  saddle  point  method. 

1.   INTRODUCTION 
Hsu  [1]  has  given  the  asymptotic  expansion 

S(N,K)  "~   (^K2)N_K[1  +  £  K~Sf  (N-K)  +  0(K_t"1)]/(N-K)!      (1) 

s=l    S 

for  Stirling  numbers  S(N,K)   of  the  second  kind,  where  f  (N-K)   are 

8 

polynomials  of  argument  N-K  and  f  (0)  =  0.   The  expansion  (1)  is 
useful  only  for  N-K  small,  as  is  indicated  in  Section  3.   We  obtain 
a  complete  asymptotic  expansion  of   S(N,K)   in  powers  of   (N+l) 
valid  for  all  K,   using  the  saddle  point  method  previously  employed 
by  Moser  and  Wyman  [2]  and  de  Bruijn  [3]  for  the  asymptotic  repre- 
sentation of  the  related  Bell  numbers. 

2.   ASYMPTOTICS  OF  S(N,K) 
A  generating  function  for   S(N,K)  is 

,eZ-l.K    f  K!  _/M   .   N-K  m 

Z       N=K  W" 


*This  work  was  partially  supported  by  the  Office  of  Naval  Research 
under  Contract  number  NR  042-286  at  the  Naval  Postgraduate  School. 
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Hence,  the  Cauchy  integral  formula  gives 

S(N'K)  =  27iKT  Jc  (eZ-l)Kz"N-1dz  (3) 

where  the  contour  C  encloses  the  origin.   Equating  the  derivative 
of  the  integrand  to  zero  gives  the  equation 

(t-Z)eZ-t  =  te_t  ,  (4) 

where   t  =  (N+l)/K,   for  the  location  of  the  saddle  point  of  the 
modulus  of  the  integrand.   The  principal  saddle  point   z  =  u   is  on 
the  positive  real  axis  with  t  -  1  <  u  <  t.   There  are  other  sub- 
sidiary saddle  points  at  complex  roots  of  (4) ,  which  we  neglect  in 
comparison  with  the  higher  saddle  point  at  z  =  u.   Since  there  are 
no  other  roots  of  (4)  for   |t-z|_^t-u,  we  may  apply  the 
Lagrange  inversion  formula  to  obtain 

u  =  t  -  t    mm-1(te-t)m/ml  (5) 

m=l 

convergent  for   t  >  1.   Another  form  of  (4)  is  the  identity 

K  =  (N+l)(l-e"U)/u  (6) 

needed  later.   Since  the  axis  of  the  saddle  point  is  perpendicular 
to  the  real  axis,  the  part  of  the  contour  C  descending  from  z  =  u 
is  taken  as  the  line  z  =  u+iy,  -<»  <  y  <  <»,  parallel  to  the 
imaginary  axis.   This  path  has  the  property  that  the  modulus  of  the 
integrand  in  (3)  is  maximal  at  the  saddle  point   z  =  u.   The  closed 
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contour  C  is  completed  by  a  half  circle  of  infinite  radius 
enclosing  the  origin.   The  contribution  to  the  integral  (3)  on  this 
semicircular  path  is  zero  since  N  >  0.   The  integral  in  (3)  then 
becomes 

/oo 
exp  ^(u+iy)dy  (7) 

—00 

where 

i|/(z)  -  K  ln[(eZ-l)/(eU-l)]  -  (N+l)  ln(z/u)  .  (8) 

The  contribution  of  the  various  parts  of  the  z  =  u  +  iy  path 
to  the  integral  must  now  be  studied.   As   |exp  i^(z)  |  =  exp  Re  iJj(z) 
we  have  to  study 

Re  iKu+iy)  -  K  ln[  (e2u-2eUcos  y+l)'1/  (eU-l)  ]-(N+l)ln(l+y2u~2)^  . 


We  shall  show  that  we  can  restrict  ourselves  essentially  to  the 
interval   |y|  <  tt.   Since 


2  -2  -2 

1  +  y  u   >_  1  +  Tr(2y—rr)u 


for  y  >  tt  we  have 


,  u  1NK  -N-l 
(e  -1)  u 


f 

J    TT 


exp  ijj(u+iy)  dy 


<  u    (e  +1)  /tt(N-I)  (1+tt  u   ) 


which  is  of  0(N  NeN)   for  K  small  and  of  0(N  e      )   for  K 
large.   Since  i^Re  i>(u+iy)   is  even  the  part  of  the  integral  (7)  for 
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|y|  >  tt  tends  toward  zero  as  N  ->  °°.   We  now  direct  our  attention 
to  the  interval   |y|  <  tt  where  the  saddle  point  at  y  =  0  gives 
the  main  contribution.   The  Taylor  expansion  of  i|;(u+iy),  convergent 
for   |y|  <  u,  is 


N+M  1      x    2        ......     f     (iy)3+2 ,  dNj+lr    1-e  U  1, 

♦  =  "  1^(u  "  TT)y     +   (N+1)    L  ik)!^        l~TT:  ~  z]z=u      (9) 
e  -1  j=l       J  u(e  -1) 


where  the  identity  (6)  has  been  used.   We  now  make  the  substitutions 


w  =  [-=-]  [1 — ]  y/u  (10) 

eU-l 


and 


a,  -  (iWu)^+2(-^)3+1[-i^-  -  ±1/0+2)!  ft  -  -^-]^+1    (11) 
j  dz      ,  z  n .    z  z=u  u  n 

u(e  -1)  e  -1 


to  obtain 


S(N,K)    .   B    fVw2+f[(N+1)~\w  (12) 


where 


B  =  N!(eU-l)K/Tr/2(N+l)K!uN/l+u/(l-exp  u)  (13) 


and 


i  i 

f[(8+i)~H]   =  Y,    a.(N+l)'^j  .  (14) 

j=l  J 

On  using  Cauchy's  theorem  for  derivatives  we  find  that 

|  a.  |  <  L  aj  (15) 

J 
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where 


and 


L  =  2w2u2/(u-l)(l-e1"u)(l  -  -^-)  (16) 

e  -1 


a  =  /2wu/(l — )"*  .  (17) 

e  -1 

-h 
We  now  expand  exp  f[(N+l)   ]   in  a  Taylor  series  of  the  form 


J*     J_         -\-  -i 

21   =        K   CMX1  "\   2J  (18) 


exp  f[(N+l)  *]   =  ;   b.(N+l) 
j -0  J 


where     b.   =   1     and     b .      are  even  or  odd   polynomials   in     w     according 
to     j      whether   it   is  even  or  odd.      By  a   lemma  of  Moser  and  Wyman    [A] 

|b.|     :  LcrUl+L^"1   .  (19) 

By  (15)  and  (19)  we  see  that  the  series  (14)  and  (18)  are  convergent 
for  suitably  large  N.   At  this  point  all  of  the  reasoning  given  by 
Moser  and  Wyman  [2]  for  the  asymptotic  expansion  of  the  Bell  numbers 
applies  equally  well  to  the   S(N,K)   expansion  and  we  conclude  that 


s-1       -'     f°°  2 

S(N,K)^B{^   (N+l)  J  i    b   e'W  dw  +  0[(N+1)'S]}  .       (20) 
j=0        J  -=°   J 


The  first  two  terms  of  the  expansion  (20)  have  been  calculated  to  be 
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S(N,K)~ 


N!(eU-l)K 


/2tt(N+1)K!u Vl+u/(l-exp   u) 


-      3,   2u.,    u,,  N  .  u   ..,-3        cro      2,  u,lW  u  .v-2,2 
;i        6-u    (e     +4e  +1)  (e   -1)        _   5[2-u    (e  4-1)  (e   -1)      ]      ■, 

8(N+l)[l+u/(l-exp  u)]2  24(N+l)[l+u/(l-exp  u)]3 


(21) 
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3.   NUMERICAL  EXAMPLE 

The  following  6-signif icant-f igure  table  compares  the  exact 
values  of   S(100,K)   with  values  computed  from  (20)  and  (1)  for 
k;   1  <  k  <  100: 


S(100, 

K) 

S(100,K) 

S(100,K) 

S(100,K) 

K 

Exact 

1  term  c 
1.00083 

»f  (20) 

2  terms  c 
1.00000 

>f  (20) 

4  terms 
7.74429 

of  (1) 

1 

1.00000 

io"176 

2 

6.33825 

io29 

6.34348 

io29 

6.33825 

io29 

1.81186 

io"115 

3 

8.58963 

io46 

8.59672 

io46 

8.58962 

io46 

3.59066 

io"80 

4 

6.69558 

io58 

6.70110 

io58 

6.69557 

io58 

2.98751 

io"55 

5 

6.57384 

io67 

6.57927 

io67 

6.57384 

io67 

4.46215 

io"36 

6 

9.07387 

io74 

9.08136 

io74 

9.07387 

io74 

1.41738 

io"20 

7 

6.41760 

io80 

6.42291 

io80 

6.41760 

io80 

1.23857 

io"7 

8 

5.05211 

io85 

5.05630 

io85 

5.05210 

io85 

1.38293 

io4 

9 

7.31910 

io89 

7.32522 

io89 

7.31910 

io89 

5.33906 

io13 

10 

2.75500 

io93 

2.75733 

io93 

2.75500 

io93 

1.41107 

io22 

11 

3.44958 

io96 

3.45255 

io96 

3.44958 

io96 

4.16533 

io29 

12 

1.72552 

io" 

1.72704 

io" 

1.72552 

io" 

1.97468 

io36 

13 

3.96430 

io101 

3.96785 

io101 

3.96430 

io101 

1.98352 

io42 

14 

4.66349 

io103 

4.66773 

io103 

4.66350 

io103 

5.24140 

io47 

15 

3.06221 

io105 

3.06502 

io105 

3.06221 

io105 

4.32888 

io52 

16 

1.20332 

io107 

1.20443 

io107 

1.20333 

io107 

1.28480 

io57 

17 

2.99575 

io108 

2.99850 

io108 

2.99576 

io108 

1.53680 

io61 

18 

4.95343 

io109 

4.95793 

io109 

4.95345 

io109 

8.15008 

io64 

19 

5.65920 

io110 

5.66423 

io110 

5.65922 

io110 

2.07645 

io68 

20 

4.61926 

io111 

4.62326 

io111 

4.61927 

io111 

2.72088 

io71 
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K 

S(100,K) 
Exact 

1 

2 

S(100,K) 
term  of  (20) 

112 
.77402  10 

S(100,K) 
2  terms  of  (20) 

112 
2.77171  10 

4 
1, 

S(10( 
terms 

.94406 

),K) 

of  (1) 

21 

2.77170 

io112 

io74 

22 

1.25295 

io113 

1 

.25396 

io113 

1.25295 

io113 

7, 

.96614 

io76 

23 

4.35876 

io113 

4, 

.36209 

io113 

4.35877 

io113 

1 

.95614 

io79 

24 

1.18873 

io114 

1, 

.18959 

io114 

1.18873 

io114 

2 

.99138 

io81 

25 

2.58320 

io114 

2, 

,58496 

io114 

2.58321 

io114 

2 

.94696 

io83 

26 

4.53753 

io114 

4, 

.54042 

io114 

4.53753 

io114 

1, 

.92723 

io85 

27 

6.52519 

io114 

6, 

,52905 

io114 

6.52520 

io114 

8, 

.59300 

io86 

28 

7.76973 

io114 

7, 

,77396 

io114 

7.76973 

io114 

2, 

.67528 

io88 

29 

7.73865 

io114 

7. 

,74251 

io114 

7.73864 

io114 

5, 

.94163 

io89 

30 

6.50629 

io114 

6. 

,50922 

io114 

6.50628 

io114 

9 

.59702 

io90 

31 

4.65568 

io114 

4, 

,65757 

io114 

4.65567 

io114 

1, 

.14719 

io92 

32 

2.85660 

io114 

2. 

,85763 

io114 

2.85659 

io114 

1, 

.03105 

93 
10 

33 

1.51311 

io114 

1, 

,51358 

io114 

1.51310 

io114 

7, 

.06841 

93 
10 

34 

6.96173 

io113 

6. 

,96362 

io113 

6.96170 

io113 

3, 

.74513 

io94 

35 

2.79791 

io113 

2, 

,79856 

io113 

2.79790 

io113 

1, 

.55216 

io95 

36 

9.87317 

io112 

9. 

,87504 

io112 

9.87311 

io112 

5, 

.08767 

IO95 

37 

3.07351 

io112 

3. 

,07397 

io112 

3.07349 

io112 

1, 

,33233 

IO96 

38 

8.47721 

io111 

8, 

,47817 

io111 

8.47715 

io111 

2, 

,81367 

IO96 

39 

2.07993 

io111 

2, 

,08009 

io111 

2.07991 

io111 

4, 

.83323 

io96 

40 

4.55640 

io110 

4, 

,55659 

io110 

4.55636 

io110 

6, 

,80712 

io96 

41 

8.94238 

io109 

8. 

,94247 

io109 

8.94230 

io109 

7, 

,91868 

io96 

42 

1.57729 

io109 

1, 

,57725 

io109 

1.57727 

io109 

7, 

,66078 

io96 

43 

2.50762 

io108 

2. 

,50750 

io108 

2.50759 

io108 

6. 

,20269 

1096 

44 

3.60314 

io107 

3. 

,60287 

io107 

3.60310 

io107 

4. 

,22802 

IO96 

45 

4.69092 

io106 

4. 

,69045 

io106 

4.69087 

io106 

2. 

,43965 

IO96 
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K 

S(100,K) 
Exact 

S(100,K) 
1  term  of  (20) 

5.54563  IO105 

S(100,K) 
2  terms  of  (20) 

5.54625  IO105 

S(10( 

4  terms 

1.19777 

),K) 
of  (1) 

46 

5.54631 

10105 

io96 

47 

5.96839 

101M 

5.96753 

io104 

5.96832 

io104 

5.02740 

io95 

48 

5.85718 

10103 

5.85623 

io103 

5.85711 

io103 

1.81206 

io95 

49 

5.25178 

io102 

5.25084 

io102 

5.25171 

io102 

5.63198 

io94 

50 

4.30983 

io101 

4.30900 

io101 

4.30977 

io101 

1.51529 

94 
10 

51 

3.24224 

10ioo 

3.24157 

io100 

3.24219 

io100 

3.54197 

93 

10" 

52 

2.23926 

io" 

2.23877 

io" 

2.23922 

io" 

7.21740 

92 
10 

53 

1.42177 

IO98 

1.42146 

io98 

1.42175 

IO98 

1.28610 

io92 

54 

8.30948 

io96 

8.30758 

io96 

8.30936 

io96 

2.01003 

io91 

55 

4.47546 

IO95 

4.47442 

io95 

4.47539 

io95 

2.76288 

1090 

56 

2.22373 

94 

io*q 

2.22322 

94 
10 

2.22370 

94 

ioy4 

3.34862 

IO89 

57 

1.02031 

io93 

1.02007 

io93 

1.02029 

io93 

3.58715 

io88 

58 

4.32678 

io91 

4.32581 

io91 

4.32671 

io91 

3.40392 

IO87 

59 

1.69718 

io90 

1.69682 

io90 

1.69716 

io90 

2.86713 

io86 

60 

6.16213 

io88 

6.16086 

io88 

6.16202 

io88 

2.14775 

IO85 

61 

2.07225 

io87 

2.07185 

io87 

2.07221 

io87 

1.43335 

IO84 

62 

6.45799 

io85 

6.45686 

io85 

6.45787 

io85 

8.53612 

io82 

63 

1.86594 

io84 

1.86566 

io84 

1.86591 

io84 

4.54310 

io81 

64 

5.00048 

io82 

4.99983 

io82 

5.00038 

io82 

2.16381 

IO80 

65 

1.24327 

io81 

1.24315 

io81 

1.24325 

io81 

9.23409 

io78 

66 

2.86851 

io79 

2.86832 

io79 

2.86846 

io79 

3.53479 

io77 

67 

6.14247 

io77 

6.14228 

io77 

6.14235 

io77 

1.21495 

IO76 

68 

1.22080 

io76 

1.22082 

io76 

1.22078 

io76 

3.75283 

io74 

69 

2.25191 

io74 

2.25204 

io74 

2.25186 

io74 

1.04256 

io73 

70 

3.85478 

io72 

3.85521 

io72 

3.85470 

io72 

2.60661 

io71 
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_K_ 

S(100,K) 

Exact 

S(100,K) 
1  term  of  (20) 

6.12310  IO70 

S(100,K) 
2  terms  of  (20) 

6.12193  107° 

S(10C 

4  terms 

5.86843 

>,K) 
of  (1) 

71 

6.12206 

1070 

io69 

72 

9.01786 

io68 

9.01999 

io68 

9.01766 

io68 

1.19026 

io68 

73 

1.23149 

1067 

1.23187 

io67 

1.23147 

io67 

2.17566 

io66 

74 

1.55828 

io65 

1.55889 

io65 

1.55825 

IO65 

3.58501 

io64 

75 

1.82584 

io" 

1.82671 

io63 

1.82579 

io63 

5.32626 

io62 

76 

1.97939 

io61 

1.98054 

io61 

1.97934 

io61 

7.13564 

io60 

77 

1.98357 

io59 

1.98495 

io59 

1.98353 

io59 

8.62061 

io58 

78 

1.83542 

io57 

1.83692 

io57 

1.83538 

io57 

9.39120 

io56 

79 

1.56619 

1055 

1.56769 

io55 

1.56615 

io55 

9.22442 

io54 

80 

1.23065 

io53 

1.23202 

io53 

1.23062 

io53 

8.16802 

io52 

81 

8.88952 

io30 

8.90092 

io50 

8.88930 

IO30 

6.51851 

io50 

82 

5.89155 

IO48 

5.90025 

io48 

5.89140 

io48 

4.68689 

io48 

83 

3.57456 

io46 

3.58062 

io46 

3.57447 

io46 

3.03474 

io46 

84 

1.98038 

io44 

1.98423 

io44 

1.98032 

io44 

1.76835 

io44 

85 

9.98901 

io41 

1.00113 

io42 

9.98874 

io41 

9.26435 

io41 

86 

4.57157 

io39 

4.58323 

io39 

4.57144 

IO39 

4.35788 

io39 

87 

1.89083 

io37 

1.89636 

io37 

1.89078 

IO37 

1.83701 

io37 

88 

7.03510 

io34 

7.05874 

io34 

7.03487 

IO34 

6.92073 

io34 

89 

2.34174 

32 
10 

2.35081 

32 

10 

2.34166 

32 

10 

2.32155 

32 
10 

90 

6.92830 

io29 

6.95936 

io29 

6.92803 

io29 
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